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In order to meet the requirements for precision pointing of 34-m antennas, ade- 
quate control design and simulation software have to be developed, with a detailed 
description of the supporting analytical tools. This article describes a control sys- 
tem model for the elevation drive of the DSS 13 antenna. The model allows one 
to simulate elevation dynamics, cross-coupled dynamics in azimuth and elevation, 
and RF pointing error. A modal state-space model of the antenna structure was 
obtained from its finite-element model with a free rotating tipping structure. Model 
reduction techniques were applied separately for the antenna model and rate-loop 
model, thereby reducing the system order to one-third of the original one while 
preserving its dynamic properties. Extensive simulation results illustrate properties 
of the model. 


I. Introduction 

In order to meet the requirements for the 34-m antenna 
pointing accuracy, an appropriate control system model 
has to be developed. The purpose of this article is to 
give a detailed description of the antenna control system 
modeling. This analytical approach is the background for 
the computer software development, which is used for the 
simulations of antenna dynamics due to control input, dis- 
turbances, and parameter variations, and is a tool for the 
control system design. 

A state-space control system model for elevation point- 
ing of the DSS 13 antenna is described. Special attention 
is paid to antenna structure modeling. A finite-element 
model of the antenna structure with a free rotating tip- 
ping structure is used to obtain its state-space model in 
modal coordinates. The model was developed with the 
following assumptions: 

(1) Linearity of the model (nonlinear effects due to dry 
friction are omitted). 


(2) The structural model consists of mode shapes up to 
5 Hz. 

(3) The elevation angle of the tipping structure is 90 deg. 

(4) The azimuth rotation joint is locked. 

(5) No disturbances act on the structure. 

The balanced reduction technique was applied to re- 
duce the structural model and subsequently the rate-loop 
model. The order of the resulting reduced model is one- 
third of the original one, while the modeling error is in- 
significantly small. Extensive simulations were performed 
and are described below to illustrate the control system 
properties. 1 


1 W. Gawronski, “Software Package for Modeling and Simulation of 
the Elevation Control System for the DSS 13 Antenna,” JPL In- 
teroffice Memorandum 3324-91-016 (internal document), Jet Pro- 
pulsion Laboratory, Pasadena, California. February 6, 1991. 
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II. Structural Model 

Here and throughout the article, the state-space repre- 
sentation of a system is used. This representation is de- 
noted by a triple ( A,B,C ), which represents a first-order 
differential equation in the form 


x - Ax + Bu, y = Cx (1) 

Assuming the state vector x of dimension nxl, the input 
u of dimension p X 1, and the output y of dimension q x 1, 
the matrices A, B, and C are of dimensions n x n, n x p, 
and q x n, respectively. 

The state-space model for an antenna structure with 
a 90-deg elevation angle and locked azimuth rotation is 
obtained. A balanced reduction technique applied to the 
structural model gives the smallest acceptable model. 

A. Full Model 

In this section the state-space triple (A , , B a ,C a ) for the 
antenna structure is determined. The triple is obtained in 
modal coordinates. The state vector of the structure x a 
consists of modal displacements q, and modal velocities 

u,: 


= [« 7 \ vj] ( 2 ) 

The antenna structural model is generated from the 
natural frequencies u>,- and mode shapes </>,-, obtained from 
the finite-element model generated by JPL/IDEAS, and 
presented in Appendix A. Modal damping of the struc- 
ture C i is assumed to be 0.5 percent, i.e., Q = 0.005 for 
i = 1 , . . . , 10. The antenna is a free rotating structure 
with a joint between the alidade and tipping structure. 
Hence, the open-loop structural model is unstable. From 
the finite-element model, m = 10 modes are obtained; 
the first one is the rigid-body mode with zero natural fre- 
quency. The modes are determined for selected p points 
of interest. Thus, <t>i = [<t>n , <f>i 2 , ■ ■ ■ , <i>ip\ T for i = 1, . . .,10. 


where I m is the identity matrix of dimensions m x m. De- 
tails of the derivation of A a are presented in Appendix B. 
The numerical values of f2, Z, and other parameters are 
given in Appendix A. 

2. Matrix B a . The bull-gear-pinion connection is 
shown in Fig. 1. The rollers of the connection are redun- 
dant; they compensate for gravity forces. The rigid pinion 
housing includes two drive systems. A torque applied at 
the pinion of the gearbox is the input to the structure. The 
pinion radius is r p = 6 inches. A node labeled no = 86302 
in the finite-element model is a point of contact between 
the pinion and the bull gear, while a node with the la- 
bel nu = 86881 is located at the joint of the supporting 
truss structure and the pinion housing. The joint nu is at 
distance p from node no. The distance p is small in com- 
parison to the bull gear radius; thus, in further analysis 
the nodes nu and no coincide. The torque T applied to 
the pinion is in equilibrium with tangential (to the bull 
gear and the pinion) forces F no and F nu , see Fig. 2. For 
the force F no with y and z components F noy and F noz , and 
the absolute value |F no |, one obtains from Fig. 2 

\F no \ = \T\/r p , 

Fnoy ~ |-^no| sin Or, F noz — l-^nol COSO 


Hence, 


Fnoy 


r p 1 sin a 

F 

_ 1 noz _ 


r ~ 1 cos a 


T = B no T 


where B n0 = r^^sina cosa] T and o = 28 deg. From 
Fig. 2, one can see that F nu = —F no ; hence, F nu = 
—B no T. In global coordinates, the force F a acting on 
the structure is a vector Fj = [0 0 • • • F% 0 0 • • • F^ u 
0 ••• 0]. Thus, 


1. Matrix A a . Denote 


F. = B 0 T 


n = diag{wi), Z = dta ff (Ci), i=l,2,...,m (3) where Bj = [0 0 - • • B% 0 0 - B 0 ••• 0]. The 

matrix B a is transformed into modal coordinates: 


Then, the system matrix A, for the antenna structure is 


A, 


0 Im 
-D 2 -2Zf2 


( 4 ) 


B, = M~ l * T B 0 (5a) 

where M m = diag{m m i, . . . , m mm ) is the modal mass ma- 
trix. Denoting <j> no = lnoi • • • I 4*nu = [^lnu)---) 
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4 >mnu ] i $ Jno — inoy <Ainoz]i an d 0inu — \ 4 >inuy <t>inuz\, 

where <t>i„ oy is the «th mode component at node no, at 
y direction (mode shapes are listed in Appendix A), the 
above equation reduces to 

B, = M- 1 ($ n0 -$ nu )5„ 0 (5b) 

3. Matrix C a . The following outputs of the antenna 
structure are considered: elevation angle 8 e , elevation rate 
8 e , pinion rate 9, pointing error angle in the x direction 
E x , and pointing error angle in the y direction e y . Hence, 
the antenna output matrix C , consists of five rows: C,\ 
for the elevation angle, C, 2 for the elevation rate, C a 3 for 
the pinion rate, C, 4 for pointing error in the x direction, 
and C a 5 for pointing error in the y direction: 

Cj = [ ci ; cj 2 ci ; cj 4 cj 5 ] ( 6 ) 

The first two rows of C, are determined as follows. In 
the finite-element model, the node at the bull gear center 
has a label nb = 5380, and the node at distance R to nb 
has a label nc = 41212 (see Fig. 2). High stiffness of the 
bull gear and the close location of the two nodes allows one 
to determine the elevation angle as a rigid-body rotation: 

9 e — (f/nc 2/nfc )/ R 

= (Cncyty C n by < &)x a / R — ( 4 >ncy < t > nby) x > / ^ 

where C ncy and C n b y denote rows with all but one zero 
elements. The nonzero element of C n cy is equal to 1, and 
is located in the position corresponding to the location of 
y displacement at node nc in the vector x a . The nonzero 
element of C„ is equal to 1, and located in the position 
corresponding to the location of y displacement at node nb 

in the vector x a . Furthermore, <j> nC y = [0incy <t>mncy ], 

where <pi ncy is the ith mode component at node nc at y 
direction, and <j>nby — \<f>inby , • • ■ , <f>mnby]i where <frinby is 
the z'th mode component at node nb at y direction. From 
the above, it follows that 

C,1 = [0 C ae ], C a 2 = [C ae 0] (7a) 

where C at — ( <t>ncy finby) / R- 

For the determination of the pinion rate measurement 
matrix C s 3, denote the velocity at pinion housing v p (at 
node nu), bull gear velocity tq at node no (see Fig. 2), 


and their y and z components by v py , v pa , Vb y , and vi u , 
respectively. Projections of uj and v p onto directions tan- 
gential to the bull gear at node no are denoted tq and tq, 
respectively. 

tq = Vby sin a + cos a, vi = v py sin a + v pz cos a 
The pinion rate is 0 = (tq — )/r p , and noting that 

Vby = </>noyV a , Vbz — t t>nozV a , V py — <t>nuyV a , and V pz — 

<t>nuiV a , where v a is a modal velocity defined in Eq. (2), one 
obtains 

tq V2 = [(^no? (frnuz ) COS Ot -f- ( <j>noy ^nuy ) sin Qrj V s 

Therefore, 9 = C ap v a = [0 C ap ]x a , giving finally 

C a 3 = [0 C ap ] (7b) 

where C ap = [( <f>noz ^nuz ) COS Ot -1- (0noy 4*nuy ) sin Of] /v p . 

The last two rows of C a , the output matrices for the 
pointing error in the x and y directions, are: 

c a 4 = [C ex 0], C a5 = [C ty 0] (7c) 

where C tx = [333.26, 0.07973, -0.60701, 0.15431, 48.635, 
4.0901, -40.343, 113.24, 225.17, -39.59] and C ey = 
[-9.0557, -24.706, -6.6233, 3.6844, -19.179, 23.571, 
43.871, 180.87, —59.909, —227.56] are obtained from the 
JPL/IDEAS. 

B. Model Reduction 

The structural model under consideration consists of 
10 modes or 20 states. The set of modes may consist of 
modes which are not important in system dynamics and 
can thus be eliminated. Observability and controllability 
properties are used to determine which modes to eliminate. 
Observability is a property of a system that indicates the 
participation of each system state in the output with all 
states equally excited. Controllability, on the other hand, 
indicates the excitation of each state due to white noise, 
or impulse input. A balanced representation [1] is a state- 
space representation with its states equally controllable 
and observable, and the Hankel singular value is the mea- 
sure of the joint controllability and observability of each 
balanced state variable. The states with small Hankel sin- 
gular values can be deleted, since they are weakly excited 
and weakly observed at the same time, causing minimal 
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modeling error. For flexible structures, modal represen- 
tation is almost balanced, assuming small damping and 
distinct poles [2,3]. Each mode has almost the same con- 
trollability and observability property, hence, each mode 
can be considered for reduction separately. The Hankel 
singular value for the *th mode is given by Eq. (53) of [2] 
and Eq. (14) of [4]: 


(b.ibji)(cj i c, i ) 


where 6 S1 - is the ith row of B, and c,i is the ith col- 
umn of C,. For the rigid-body mode, the Hankel singular 
value does not exist — its value tends to infinity — hence, 
the rigid-body mode is always included in the reduced 
model. For each mode, the Hankel singular value is de- 
termined to help the designer decide on the number of 
modes in the reduced structural model. 

Hankel singular values of the 10-mode model of the an- 
tenna ( A,,B,,C ,) are plotted in Fig. 3, where the Hankel 
singular value of the rigid-body mode is omitted. The 
reduced-order model consists of five modes: the rigid- 
body mode and four elastic modes with natural frequen- 
cies 3.1358, 4.1149, 4.4503, and 4.5688 Hz, indicated in 
Fig. 3. The transfer functions of the full and reduced mod- 
els of the antenna are shown in Fig. 4, indicating that the 
reduced model fully represents properties of the original 
system. Note that disturbances are not included in the 
model; therefore, they are not considered in model reduc- 
tion. One should note, however, that a different reduced 
model would be obtained for a model with disturbances 
included. 


III. Rate-Loop Model 

This section presents details of rate-loop modeling and 
introduces model reduction. 

A. Full Model 

The rate-loop model of the antenna is shown in Fig. 5. 
It consists of two identical subsystems denoted G 0 , an- 
tenna structure G a , command amplifier G c , and the bias 
loop, with the bias command Vbi a ,. 

The subsystem G c is shown in Fig. 6 (a). The first state 
variable is denoted = iq and its state equation is 

X{ = -rf 1 *! + kik, rf 1 9 0 (9) 


The subsystem G 0 , Fig. 6 (b), consists of four subsystems 
itself: G\, G 2 (both amplifiers), G 3 (motor armature), and 
G 4 (motor and gearbox mechanical properties). The state 
equations for G 1 are 


X2 = *3 

( 10 a) 

X 3 — X\ 7Vj £3 ktach% 7 

( 10 b) 

where £7 = u m . Similarly, for G 2 ’ 


£4 — £5 

( 11 a) 

Z 5 = k r T^ 1 X 2 + k r T 2 T^ l X 3 - T$ l X 5 - k c ur X 6 

(lib) 

where X 6 = i 0 . For G 3 , from Fig. 6 : 


i 6 = k f k i L~ 1 T^ 1 X 4 + k f k i T 4 L~ 1 Ts l x 5 


- R a L~ l x e - k i L~ 1 x 7 

( 12 a) 

To — k m XQ 

( 12 b) 


The system G 4 denotes the motor and gearbox model. 
There are two inputs to the system: motor torque T 0 and 
pinion rate 6 . The output is the pinion torque T. The 
gearbox flexibility is described either by the input or the 
output stiffness, as shown in Figs. 7(a) and 7(b). For the 
output stiffness model, with notations as in Fig. 7(a), one 
obtains equations for the gearbox rate 

T n = NT a , e N = 9 m /N 

and the dynamic equilibrium equations 
J m u m + Tn /N = T 0 

along with the spring equation for the output stiffness k go 

T = T n = k go (0 m - N9)/N 

Similar equations are derived for the input stiffness 
model in Fig. 7(b). Both models are represented by Fig. 8 , 
where N x = l/7Vand k gx = k go for the output stiffness 
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model and N x = N and k gx = k g i for the input stiffness 
model. Denoting x§ = T, the state equations are 

X 7 = k m Jm l X 6 - J~ 1 N~ 1 x 8 (13a) 

x 8 = k gx N x x 7 - k gx N x NC„x, (13b) 

where N x = N for k gx = k g i, and N x = \/N for k gx - k g0 . 

Combining Eqs. ( 10)— ( 13) for models G\, G 2 , G 3 , 
and G 4 , one obtains the state-space model for G 0 (in- 


puts v\ and 9, outputs T and i„, and state vector x 0 = 

[•^ 2 ) X3, X4, X5 , Xg , X7, Xg] ). 


xo -(- Bum 4- B \20 

(14a) 

T 

— C 1 x 0 

*0 

(14b) 


where 


A 1 = 


0 

0 

0 


0 

0 


1 

-I/7-3 

0 


k r /r 3 k r T 2 /T 3 

0 0 


0 

0 

1 

-I/75 


0 

0 

0 

-k 


0 

ktach 

0 

0 


kfki/L a T 5 kj ki T 4 / L a Rq/Lq kb/ L a 


km/ Jn 

0 


0 


k gx N x 


0 

0 

0 

0 

0 

-1 /JmN 

0 


sfi = [0 1 0 0 0 0 0], B '[ 2 = [0 0 0 0 0 0 - k gx N x N], C n = [0 0 0 0 0 0 1], 


C 12 = [000 0 1 0 0], Ci = 


Cn 

C 12 


B\ = [B\\ B u \ 


Two systems’ Go’s are located in the rate-loop system 
(see Fig. 5): the first with inputs Un and 0, outputs T\ 
and i 0 i, and state vector X21, and the second with inputs 
V 12 and 9, outputs T 2 and i o2 , and state vector X22- The 
equations for these systems are 


*21 — ^1*21 + •fiiit'11 + B \20 (15a) 


T 2 

io2 


= C1X22 


(16b) 


The subsystem G c is an amplifier for the bias signal, with 
input Av c and output v u . Denoting x„ = v u , its state 
space equation is obtained: 


'T: 

*01 


— C1X21 


and 


(15b) 


X a = -T 6 1 X a + kctfrTe 1 Av c (17) 


The antenna structure equations according to Section II 
are 


X22 — ^1^22 + B 11U12 + B \20 


(16a) 


x, = A s x, + B S T, 9 e = C, 2X 5 , 9 = C, 3 x, (18) 
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The nodal equations are 

T = Ti -f T 2 , Ti = Cn* 2 i> ^2 = C 11 I 22 (19a) 

Au c = V c i — V C 2 + kb s Viia, = k ct C 12(2:21 — X22) + ki,,Vi,i a> 

(19b) 

t>n = ui — t;„ = xi — x a , ^12 = tq + v u = xi + x a 

(19c) 

Denoting the state vector xj = [x^, x ^ x , x^ 2 , x a > X J] and 
combining Eqs. (9)-(19), the equations for the whole sys- 
tem can be written 

x r = A T x T + B r u, y = C T X T (20) 


balancing procedure from [3] is used, and is summarized in 
Appendix C. The procedure from Section II is applicable 
to flexible structures only. 

The rate-loop system, with already reduced antenna 
structural modes, has 26 state variables. A plot of the sys- 
tem Hankel singular values is shown in Fig. 9. The Hankel 
singular values of component number 15 and larger are 
small in comparison to the remaining components. There- 
fore, the reduced system consists of 14 state variables. The 
transfer functions of the reduced and full rate-loop systems 
are compared in Fig. 10. The full system contains the 
full structural model: 36 state variables altogether. From 
Fig. 10 one can see that the reduced model of 14 state vari- 
ables almost exactly approximates the full 36-state model. 
In Fig. 11 step responses of full and reduced models are 
compared, showing good agreement between the full and 
reduced models. 


where u T ■■ 

= [0o Uiiaj] and y T 

= [61 b T t 9 T ], and 


'-l/n 

0 


0 

0 

0 


Bn 

A\ 


0 

-B n 

BnC,3 

A t = 

B n 

0 


A\ 

Bn 

B\2C,z 


0 

Pi 


-Pi 

-l/ffc 

0 


. 0 

B,C\\ B s Cn 

0 

A, 


kik s /n 

0 

0 0 

O' 



Bj = 




» 



0 

0 

O 

to 

0 




■000 

0 

C ti - 




C r = 

0 0 0 

0 

C, 2 

J 




.0 0 0 

0 

c, 3. 





Pi = k c tfrk c tCi2/Te, P 2 = k ct f r k c t/T6 


IV. Position-Loop Model 

The rate-loop system, with the position loop closed, is 
shown in Fig. 12. A proportional-plus-integral (PI) con- 
troller is applied. For the state x 0 defined x ^ = [xj x,], 
the state-space triple of the series connection of the rate- 
loop system and the controller is 



\A r B t ku 1 


Bt kpp 

A 0 = 

1 

O 

O 

1 

, B 0 = 

1 


Co = [Ct 0 ], 


Bob = 


Bn 

0 


( 21 ) 


where k pp and kn are the proportional and integral param- 
eters of the controller, B t is the first column of B r , Bn is 
the second column of B r , and Ct is the first row of C r . 
The closed-loop system matrix is 

A c , = A 0 - B 0 C 0 (22) 


B. Model Reduction 

The rate-loop system model is reduced by applying the 
balancing principle. The system is transformed into bal- 
anced coordinates where each state is equally controllable 
and observable and not correlated with other states. The 


and the closed-loop triple is (A c i, B 0 ,C 0 ). The closed-loop 
system performance is illustrated in Fig. 13 in time domain 
and in Fig. 14 in frequency domain. One can see that the 
x component of the RF pointing error is much larger than 
the component in the y direction, and the y component of 
subreflector vertex displacement is much larger than the 
component in the x direction. In fact, for a symmetric 
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structure the RF y component and x displacement of the 
subreflector vertex should be zero. 

Finally, the action of the bias input vn a > is tested. Two 
step inputs axe introduced to the closed-loop system: a 
vti a , step of amplitude 2.5 V at the moment t — 0 sec and 
a command step 0 O of amplitude 0.1 deg at the moment 
t = 10 sec. The results are shown in Figs. 15(a) through 
15(c). For the first 10 sec, when the command 9„ is zero, 
the torques T\ and Ti are equal but with opposite sign, so 
that the resulting torque T acting on the structure is zero: 
v bias does not influence the structure. Nevertheless, both 
motors axe active, their rates w m i and w m 2 are nonzero, 
plotted in Fig. 15(c). After 10 sec, when the command is 
active, the system behavior is similar to the one without 
bias. 


V. Conclusions 

In this article, modeling techniques for the DSS 13 an- 
tenna structure and its rate-loop control system have been 


presented. Through balancing the system controllability 
and observability properties, a reduced system model has 
been obtained. Antenna dynamics, rate-loop dynamics, 
and dynamics of a system with a closed position loop have 
been simulated. The model allows one to simulate ele- 
vation dynamics, cross-coupled dynamics in azimuth and 
elevation, and a dynamic RF pointing error. 

Section I contains restrictions that need to be relaxed. 
The structural model could include modes of up to 10 Hz, 
with free elevation and azimuth joints. The elevation and 
azimuth control systems could be simulated jointly, expos- 
ing the cross-coupling of the two systems. Disturbances 
(wind and measurement noise) could be included in the 
model. The system could be modeled for several elevation 
angles of the range of [0 deg, 90 deg], and the dependence 
of the system dynamics on elevation angle investigated. 
Nonlinear effects due to dry friction and saturations could 
be included. Finally, field measurements need to be per- 
formed, and the model parameters adjusted such that a 
satisfactory agreement between simulations and field data 
is achieved. 
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Fig. 14. Step responses for the closed-loop system; position loop closed, k pp — 0.5, 
kjj = 0.1; input: elevation angle, rad. 
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Fig. 15. Step responses of the closed-loop system to the command and bias inputs 
step Vhjas — 2.5 V at t = 0 and step 9 0 = 0.1 deg at t = 10 sec. 






Appendix A 

Parameters of the DSS 13 Antenna 


I. Structural Parameters 

The structural parameters for the DSS 13 antenna are: 
q = 28 deg 

Ci = 0.005, i = 1, .... 10 

Q = diag( 0, 12.4482, 12.5833, 13.3612, 19.7028, 
20.1747, 25.2584, 25.8547, 27.9621, 28.7060) 

Mm = dia(/(215.4239, 223.8522, 38.1610, 17.1415, 
243.8451, 58.5778, 95.1770, 804.4224, 18.4959, 
18.3158) 

Node numbers: nc = 5380, no = 86302, nu = 86881, 
and nb = 41212 

Modal matrix: $ = [<pi, (p 2 , . . • , <Pio] 

(pi = ith mode shape, <pi 

= [<A» 1 1 <Pi2, •••> <Pii2] T 


where: 

<pn — ith mode component, node nc,x direction 
<pi 2 = ith mode component, node nc, y direction 
(pig = ith mode component, node nc,z direction 
<t>i 4 = ith mode component, node no, x direction 
<j>i 5 = ith mode component, node no, y direction 
(f>iQ = ith mode component, node no, z direction 
<t>i 7 = ith mode component, node nu,x direction 
<pig = ith mode component, node nu, y direction 
<pi 9 = ith mode component, node nu, z direction 
<pi io = ith mode component, node nb, x direction 
<pm = ith mode component, node nb, y direction 
<Pni — ith mode component, node nb, z direction 


The modes (in MATLAB notation): 


<P T 1 = 

[0.00000012, 

-0.00000060, 

0.00000003,... 

% 

(node 

nc, £, t/, 2 components) 


0.00000000, 

0.17696733, 

0.34187825,... 

% 

(node 

no, x, y, z components) 


0.00000000, 

-0.00000005, 

-0.00000010 ,... 

% 

(node 

nu, x,y,z components) 


0.00000012, 

-0.04215388, 

0.00000003]; 

% 

(node nb, x,y,z components) 

4>2 - 

[0.16089562, 

-0.04828905, 

-0.01058891,... 

% 

(node 

nc, x,y,z components) 


0.00000000, 

-0.00320542, 

0.00104308,... 

% 

(node 

no, x,y,z components) 


0.00000000, 

0.00012084, 

0.00023455,... 

% 

(node 

nu, x, y, z components) 


0.18123627, 

-0.05316486, 

-0.01040931]; 

% 

(node 

nb, x, y, z components) 

U>-j 

11 

[-0.02516621, 

-0.07031944, 

0.00690156,.. 

% 

(node 

nc, x,y,z components) 


0.00000000, 

-0.00103347, 

0.00003283,.. 

% 

(node 

no, x,y,z components) 


0.00000000, 

-0.00024008, 

-0.00046601,.. 

% 

(node 

nu, x, y, z components) 


-0.02842479, 

-0.07567856, 

0.00688137]; 

% 

(node 

nb, x,y,z components) 

II 

[0.00351127, 

0.03635180, 

-0.00294988,.. 

% 

(node 

nc, x, y, z components) 


0.00000000, 

0.00089782, 

-0.00009626,.. 

% 

(node 

no, x, y, z components) 


0.00000000, 

0.00012046, 

0.00023381,.. 

% 

(node 

nu, x, y, z components) 


0.00397642, 

0.03904565, 

-0.00295351]; 

% 

(node 

nb, x, y, z components) 

cn 

II 

[0.01726125, 

0.21422039, 

-0.02128878,.. 

% 

(node 

nc, x,y, z components) 


0.00000000, 

0.41745298, 

-0.03668607,.. 

% 

(node 

no, x, y, z components) 


0.00000000, 

0.01317732, 

0.02557758,.. 

% 

(node 

nu, x, y, z components) 


0.01343854, 

0.25263344, 

-0.02139128]; 

% (node 

nb, x, y, z components) 
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II 

E*ko 

[0.08392051, 

0.00028101, 

-0.03029582,... 

% 

(node 

nc, x,y,z components) 


0.00000000, 

-0.00560731, 

-0.00120284,. . . 

% 

(node 

no, x, y, z components) 


0.00000000, 

-0.00020205, 

-0.00039219,... 

% 

(node 

nu, x, y, z components) 


0.08607231, 

-0.00154962, 

-0.03083227]; 

% 

(node 

nb, x, y , z components) 

4*7 = 

[-0.02218441, 

0.16344080, 

-0.01364561,... 

% 

(node 

nc, x, y, z components) 


0.00000000, 

-0.02133899, 

0.00544590,... 

% 

(node 

no, x,y,z components) 


0.00000000, 

0.00065982, 

0.00128072,... 

% 

(node 

nu, x, y, z components) 


[-0.01743287, 

0.15653900, 

-0.01357982]; 

% 

(node 

nb, x, y, z components) 

001-J 

II 

[-0.30854177, 

0.04674442, 

0.18580163,... 

% 

(node 

nc, x,y,z components) 


0.00000000, 

0.14007292, 

0.79574209,. . . 

% 

(node 

no, x, y, z components) 


0.00000000, 

0.05702746, 

0.11069132,... 

% 

(node 

nu, x, y, z components) 


[-0.23322315, 

0.05969981, 

0.18902462]; 

% 

(node 

nb, x, y, z components) 

II 

hia> 

[-0.00970566, 

-0.01171182, 

0.00170980,... 

% 

(node 

nc, x,y, z components) 


0.00000000, 

-0.02275614, 

0.02106441,... 

% 

(node 

no, x, y, z components) 


0.00000000, 

-0.00178816, 

-0.00347085,. . . 

% 

(node 

nu, x, y, z components) 


[-0.01092381, 

-0.01958129, 

0.00171795]; 

% 

(node 

nb, x, y, z components) 

^0 — 

[-0.05660426, 

0.01088616, 

0.01266761,... 

% 

(node 

nc, x,y,z components) 


0.00000000, 

0.00316858, 

0.00525204,... 

% 

(node 

no, x, y, z components) 


0.00000000, 

0.00090903, 

0.00176443,... 

% 

(node 

nu, x, y, z components) 


-0.05696018, 

0.00916372, 

0.01287218]. 

% 

(node 

nb, x,y,z components) 


II. Rate-Loop System Parameters 

The rate-loop system parameters for the DSS 13 an 
tenna are: 

*i = 716.197 V sec/rad 
k m = 15.72 lb/A 
kb = 1.79 V sec/rad 
k, = 0.8 V/V 

ktach = 0.0384123 V sec/rad 

k r = 80 V/sec V, range: 49 — 83 
ki = 87.13 V/sec V 
ku = 0.1 
k pp =0.5 

k cur = 0.12658 V/A 
kf = 54 V/V 
k go = 1.5 x 10 7 lb/rad 


kct/r = 0.33 
k ct = 0.11111 
fcfc , = 0.66 
T\ = 0.0063662 sec 
T 2 = 0.094 sec 
t 3 = 0.002 sec 
r 4 = 0.00484 sec 
r 5 = 0.0021 sec 
r 6 = 0.7304 sec 
N = 354 

J m = 1.236 lb/sec 2 
R a = 0.456 fi 
L a = 0.011 H 
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Appendix B 

From the Finite-Element Model to the State-Space Model 


The state-space model of a flexible structure is obtained 
from its finite-element model, which consists of the mass 
M (m x m), stiffness I< (m x m), input B 0 (m x s ), output 
C og (r x m), C ov (r x m) matrices, the input u(t ) (s x 1), 
and output y(t) (r x 1). The input-output relationship is 
given by the second-order differential equation 


Denote Af” 1 /^ = fi 2 , where 0 is a diagonal ( pxp ) ma- 
trix of natural frequencies (rad/sec). At this stage a damp- 
ing matrix Z is introduced, Z = diagfa), i = 1, . . . ,p, such 
that 2 ZQ = M^Dm, and D m is a modal damping matrix 
(assumed to be known), so that from Eq. (B-5c) the modal 
model is acquired: 


M'q + Kq = B 0 u, y = C oq q + C ov q (B-l) 

where q is the vector of structural displacements. 

Consider now a modal matrix $ (m x p), p < m, which 
consists of p eigenvectors 0,- (mode shapes), i = 1 , . . . ,p: 


q m +2 ZQq m +&q m = M m l <t> T B 0 u, 


V — Coq*&Q m + C ov ft>q„ 


(B-6) 


Define a state variable x 


$ = - ■■ ,<i> P \ 

which diagonalize M and K: 


= $ T M$, K m = <t> T I<<*> 


(B-2) 


(B-3) 


X = 


’*l" 


Qm 

- X 2. 


_ Qm _ 


(B-7) 


then, Eq. (B-6) can be presented as a set of first-order 
i.e., M m and K m are diagonal (p x p) matrices of modal equations: 
mass and stiffness. When a new variable q m (p x 1) is 
introduced such that 


q = A>q m (B-4) 

and left-multiplying Eq. (B-l) by $ T , one obtains either 

<t> T M$q m + <t> T K$q m = $ t B 0 u, 

y — Coq$qm d~ C ov $q, 


Xi = x 2 

x 2 = — D 2 xi — 2Zflx2 + M~ l $B a u 
y — Coq*&% l T C ov <j)X 2 
(B-5a) or i n fhe following form 

x — Ax + Bu, y = Cx 

where 


(B-8a) 


Qm "I" KmQm — ^ V — Coq*&Qm “f" Cov^Qm 

(B-5b) 


or 



0 

i 


0 

A = 



, b = 



-n 2 

-2Zfi 


_m- 1 a> t Bo_ 


(B-8b) 


q m + M m x K m q m = M~ 1 $ T 5 0 u, 


V — CoqQqm T Cov^qn 


(B-5c) 


C = [C oq d> C ov <i>) 

is the sought state-space model in modal coordinates. 
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Appendix C 

A Balanced Representation of a Linear System 


Consider a stable, linear, time-invariant system: 

x = Ax + Bu, y = Cx, x(0) = x o (C-l) 

where x 6 R n is the system state, u € R p is the sys- 
tem input, y G R m is the system output, and (A,B,C) 
is the system state-space representation. The system con- 
trollability and observability grammians W c and W„ are 
solutions of the Lyapunov equations 

AW e + W e A T + BB t = 0, A t W 0 + W 0 A + CC t = 0 

(C-2) 

The system representation is balanced if its controllabil- 
ity and observability grammians are diagonal and equal. 
Hence, for the balanced representation (A/,, Bt,Ct) = 
(T~ l AT,T~ l B,CT), the following is true: 

w e = w 0 = r 2 , 

r = diag(-n,...,y n ), 7< > o, * = l,...,n (C-3) 


where T is a linear transformation and 7 ,• is the ith Hankel 
singular value of the system. 

The transformation T is determined as follows: 

T=PUT~\ T- 1 =r~ 1 V T Q (C-4) 

The matrices F, V, and U are obtained from the singular 
value decomposition of the matrix H: 

H = Vr 2 U T (C-5) 

where H = QP, and P and Q form the following decom- 
position of grammians: 

W c = PP T , W 0 = Q t Q (C-6) 

for example, Cholesky, or singular- value decomposition. 
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